Unveiling the Oxazolidine Character of Pseudoproline Derivatives by Automated Flow Peptide Chemistry.

Pseudoproline derivatives such as Thr(ΨPro)-OH are commonly used in peptide synthesis to reduce the likelihood of peptide aggregation and to prevent aspartimide (Asi) formation during the synthesis process. In this study, we investigate notable by-products such as aspartimide formation and an imine derivative of the Thr(ΨPro) moiety observed in flow peptide chemistry synthesis. To gain insight into the formation of these unexpected by-products, we design a series of experiments. Furthermore, we demonstrate the oxazolidine character of the pseudoproline moiety and provide plausible mechanisms for the two-way ring opening of oxazolidine leading to these by-products. In addition, we present evidence that Asi formation appears to be catalyzed by the presence of the pseudoproline moiety. These observed side reactions are attributed to elevated temperature and pressure; therefore, caution is advised when using ΨPro derivatives under such harsh conditions. In addition, we propose a solution whereby thermodynamically controlled Asi formation can be kinetically prevented.

Polymorphic amyloid nanostructures of hormone peptides involved in glucose homeostasis display reversible amyloid formation.

A large group of hormones are stored as amyloid fibrils in acidic secretion vesicles before they are released into the bloodstream and readopt their functional state. Here, we identify an evolutionarily conserved hexapeptide sequence as the major aggregation-prone region (APR) of gastrointestinal peptides of the glucagon family: xFxxWL. We determine nine polymorphic crystal structures of the APR segments of glucagon-like peptides 1 and 2, and exendin and its derivatives. We follow amyloid formation by CD, FTIR, ThT assays, and AFM. We propose that the pH-dependent changes of the protonation states of glutamate/aspartate residues of APRs initiate switching between the amyloid and the folded, monomeric forms of the hormones. We find that pH sensitivity diminishes in the absence of acidic gatekeepers and amyloid formation progresses over a broad pH range. Our results highlight the dual role of short aggregation core motifs in reversible amyloid formation and receptor binding.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion.

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Inhibitor Design Strategy for Myostatin: Dynamics and Interaction Networks Define the Affinity and Release Mechanisms of the Inhibited Complexes.

Myostatin, an important negative regulator of muscle mass, is a therapeutic target for muscle atrophic disorders such as muscular dystrophy. Thus, the inhibition of myostatin presents a strategy to treat these disorders. It has long been established that the myostatin prodomain is a strong inhibitor of the mature myostatin, and the minimum peptide of the prodomain-corresponding to the α1-helix of its lasso-region-responsible for the inhibitory efficiency was defined and characterized as well. Here we show that the minimum peptide segment based on the growth differentiation factor 11 (GDF11), which we found to be more helical in its stand-alone solvated stfate than the similar segment of myostatin, is a promising new base scaffold for inhibitor design. The proposed inhibitory peptides in their solvated state and in complex with the mature myostatin were analyzed by in silico molecule modeling supplemented with the electronic circular dichroism spectroscopy measurements. We defined the Gaussian-Mahalanobis mean score to measure the fraction of dihedral angle-pairs close to the desired helical region of the Ramachandran-plot, carried out RING analysis of the peptide-protein interaction networks and characterized the internal motions of the complexes using our rigid-body segmentation protocol. We identified a variant-11m2-that is sufficiently ordered both in solvent and within the inhibitory complex, forms a high number of contacts with the binding-pocket and induces such changes in its internal dynamics that lead to a rigidified, permanently locked conformation that traps this peptide in the binding site. We also showed that the naturally evolved α1-helix has been optimized to simultaneously fulfill two very different roles: to function as a strong binder as well as a good leaving group. It forms an outstanding number of non-covalent interactions with the mature core of myostatin and maintains the most ordered conformation within the complex, while it induces independent movement of the gate-keeper ß-hairpin segment assisting the dissociation and also results in the least-ordered solvated form which provides extra stability for the dissociated state and discourages rebinding.

Acetyl group for proper protection of ß-sugar-amino acids used in SPPS.

The synthesis of D-glucosamine-1-carboxylic acid based ß-sugar amino acids (ß-SAAs) is typically performed in nine consecutive steps via an inefficient OAc → Br → CN conversion protocol with low overall yield. Here, we present the improved and more efficient synthesis of both Fmoc-GlcAPC-OH and Fmoc-GlcAPC(Ac)-OH, ß-SAAs consisting of only 4-5 synthetic steps. Their active ester and amide bond formation with glycine methyl ester (H-Gly-OMe) was completed and monitored by 1H NMR. The stability of the pyranoid OHs protecting the acetyl groups was investigated under three different Fmoc cleavage conditions and was found to be satisfactory even at high piperidine concentration (e.g. 40%). We designed a SPPS protocol using Fmoc-GlcAPC(Ac)-OH to produce model peptides Gly-ß-SAA-Gly as well as Gly-ß-SAA-ß-SAA-Gly with high coupling efficiency. The products were deacetylated using the Zemplén method, which allows the hydrophilicity of a building block and/or chimera to be fine-tuned, even after the polypeptide chain has already been synthesized.

Covalent fragment mapping of KRasG12C revealed novel chemotypes with in vivo potency.

G12C mutant KRas is considered druggable by allele-specific covalent inhibitors due to the nucleophilic character of the oncogenic mutant cysteine at position 12. Discovery of these inhibitors requires the optimization of both covalent and noncovalent interactions. Here, we report covalent fragment screening of our electrophilic fragment library of diverse non-covalent scaffolds equipped with 40 different electrophilic functionalities to identify fragments as suitable starting points targeting Cys12. Screening the library against KRasG12C using Ellman's free thiol assay, followed by protein NMR and cell viability assays, resulted in two potential inhibitor chemotypes. Characterization of these scaffolds in in vitro cellular- and in vivo xenograft models revealed them as promising starting points for covalent drug discovery programs.

[Farnesyltransferase inhibitors have antitumoral effects in mutant KRAS containing cancer cells in preclinical models]. / Farnezil-transzferáz-inhibitorok (FTI) szinergista módon potencírozzák a mutáns KRAS inhibitorai daganatellenes hatásait preklinikai modellekben.

In silico studies raised the possibility that farnesyltransferase inhibitors (FTIs) may have antitumoral effects on KRAS mutant cancer cells. Accordingly, we have tested FTIs (tipifarnib and lonafarnib) in G12C mutant human cancer cell lines in vitro and in vivo. We have discovered that the combination of the two drugs has a synergistic antitumoral effect. Next, we have tested FTIs on G12D mutant human cancer cell lines and found that the combination has antitumoral effect in various preclinical cancer models. At last, we have also tested FTIs on G12V mutant human cancer cells and again we have detected antitumoral effects. We suggest that FTIs may have clinical relevance outside the HRAS mutant cancers.

Succinct Amyloid and Nonamyloid Patterns in Hexapeptides.

Hexapeptides are widely applied as a model system for studying the amyloid-forming properties of polypeptides, including proteins. Recently, large experimental databases have become publicly available with amyloidogenic labels. Using these data sets for training and testing purposes, one may build artificial intelligence (AI)-based classifiers for predicting the amyloid state of peptides. In our previous work (Biomolecules 2021, 11, 500), we described the Support Vector Machine (SVM)-based Budapest Amyloid Predictor (https://pitgroup.org/bap). Here, we apply the Budapest Amyloid Predictor for discovering numerous amyloidogenic and nonamyloidogenic hexapeptide patterns with accuracy between 80% and 84%, as surprising and succinct novel rules for further understanding the amyloid state of peptides. For example, we have shown that for any independently mutated residue (position marked by "x"), the patterns CxFLWx, FxFLFx, or xxIVIV are predicted to be amyloidogenic, while those of PxDxxx, xxKxEx, and xxPQxx are nonamyloidogenic. We note that each amyloidogenic pattern with two x's (e.g.,CxFLWx) describes succinctly 202 = 400 hexapeptides, while the nonamyloidogenic patterns comprising four point mutations (e.g.,PxDxxx) give 204 = 160 000 hexapeptides in total. We also examine the restricted substitutions for positions "x" from subclasses of proteinogenic amino acid residues; for example, if "x" is substituted with hydrophobic amino acids, then there exist patterns containing three x's, like MxVVxx, predicted to be amyloidogenic. If we can choose for the x positions any hydrophobic amino acids, except the "structure breaker" proline, then we get amyloid patterns with five x positions, for example, xxxFxx, each corresponding to 32 768 hexapeptides. To our knowledge, no similar applications of artificial intelligence tools or succinct amyloid patterns were described before the present work.

Synthesis of the extracellular domain of GLP-1R by chemical and biotechnological approaches.

The extracellular domain of the glucagon-like peptide-1 receptor, GLP-1R, is responsible for the binding of GLP-1, and a handful of additional agonists (such as exenatide, lixisenatide, and liraglutide) used daily for treating type II diabetes mellitus. Lead discovery and optimization, however, require binding studies, which, in turn, necessitate the total synthesis of GLP-1R, comprising 108 residues. A protein domain of 10-15 kDa size could be obtained either by expression in E. coli or by ligating solid-phase peptide synthesis (SPPS)-made fragments. However, direct overexpression fails to give a properly folded protein, as GLP-1R forms an inclusion body, which fails to refold due to improper disulfide pairing. Several bacterial strains, constructs, and fusion partners were probed and it was found that only co-expression with MBP gave a 3D-fold allowing the native disulfide bond pattern formation. Some fusion partners can act as covalently linked or in situ chaperones for guiding the refolding of GLP-1R toward success. Therefore, the bottleneck to preparing GPCR extracellular domains is the correct pairing of the Cys residues. As a proof-of-concept model, nGLP1-R was made by SPPS to form the purified full-length polypeptide chain, subjected to self-guided or spontaneous Cys pairing. However, the formation of correct SS-pairs was lagging behind any protocol in use support, and the bottleneck of large-scale protein production relies on the risky step of proper refolding, which is sometimes possible only if a suitable fusion partner effectively helps and catalysis of the correct disulfide formation.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants.

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg2+ -free state characterized by shortened ß2/ß3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2.

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

Directed Evolution-Driven Increase of Structural Plasticity Is a Prerequisite for Binding the Complement Lectin Pathway Blocking MASP-Inhibitor Peptides.

MASP-1 and MASP-2 are key activator proteases of the complement lectin pathway. The first specific mannose-binding lectin-associated serine protease (MASP) inhibitors had been developed from the 14-amino-acid sunflower trypsin inhibitor (SFTI) peptide by phage display, yielding SFTI-based MASP inhibitors, SFMIs. Here, we present the crystal structure of the MASP-1/SFMI1 complex that we analyzed in comparison to other existing MASP-1/2 structures. Rigidified backbone structure has long been accepted as a structural prerequisite for peptide inhibitors of proteases. We found that a hydrophobic cluster organized around the P2 Thr residue is essential for the structural stability of wild-type SFTI. We also found that the same P2 Thr prevents binding of the rigid SFTI-like peptides to the substrate-binding cleft of both MASPs as the cleft is partially blocked by large gatekeeper enzyme loops. Directed evolution removed this obstacle by replacing the P2 Thr with a Ser, providing the SFMIs with high-degree structural plasticity, which proved to be essential for MASP inhibition. To gain more insight into the structural criteria for SFMI-based MASP-2 inhibition, we systematically modified MASP-2-specific SFMI2 by capping its two termini and by replacing its disulfide bridge with varying length thioether linkers. By doing so, we also aimed to generate a versatile scaffold that is resistant to reducing environment and has increased stability in exopeptidase-containing biological environments. We found that the reduction-resistant disulfide-substituted l-2,3-diaminopropionic acid (Dap) variant possessed near-native potency. As MASP-2 is involved in the life-threatening thrombosis in COVID-19 patients, our synthetic, selective MASP-2 inhibitors could be relevant coronavirus drug candidates.

Solution Structure and Acid-Base Properties of Reduced α-Conotoxin MI.

The reduced derivative of α-conotoxin MI, a 14 amino acid peptide is characterized by NMR-pH titrations and molecular dynamics simulations to determine the protonation constants of the nine basic moieties, including four cysteine thiolates, and the charge-dependent structural properties. The peptide conformation at various protonation states was determined. The results show that the disulfide motifs in the native globular α-conotoxin MI occur between those cysteine moieties that exhibit the most similar thiolate basicities. Since the basicity of thiolates correlates to its redox potential, this phenomenon can be explained by the higher reactivity of the two thiolates with higher basicities. The folding of the oxidized peptide is further facilitated by the loop-like structure of the reduced form, which brings the thiolate groups into sufficient proximity. The 9 group-specific protonation constants and the related, charge-dependent, species-specific peptide structures are presented.

Synthesis of chimera oligopeptide including furanoid ß-sugar amino acid derivatives with free OHs: mild but successful removal of the 1,2-O-isopropylidene from the building block.

Complementary to hydrophobic five membered ring ß-amino acids (e.g. ACPC), ß-sugar amino acids (ß-SAAs) have found increasing application as hydrophilic building blocks of foldamers and α/ß chimeric peptides. Fmoc-protected ß-SAAs [e.g. Fmoc-RibAFU(ip)-OH] are indeed useful Lego elements, ready to use for SPPS. The removal of 1,2-OH isopropylidene protecting group increasing the hydrophilicity of such SAA is presented here. We first used N3-RibAFU(ip)-OH model compound to optimize mild deprotection conditions. The formation of the 1,2-OH free product N3-RibAFU-OH and its methyl glycoside methyl ester, N3-RibAFU(Me)-OMe were monitored by RP-HPLC and found that either 50% TFA or 8 eqv. Amberlite IR-120 H+ resin in MeOH are optimal reagents for the effective deprotection. These conditions were then successfully applied for the synthesis of chimeric oligopeptide: -GG-X-GG- [X=RibAFU(ip)]. We found the established conditions to be effective and-at the same time-sufficiently mild to remove 1,2-O-isopropylidene protection and thus, it is proposed to be used in the synthesis of oligo- and polypeptides of complex sequence combination.

Assessment of Tractable Cysteines for Covalent Targeting by Screening Covalent Fragments.

Targeted covalent inhibition and the use of irreversible chemical probes are important strategies in chemical biology and drug discovery. To date, the availability and reactivity of cysteine residues amenable for covalent targeting have been evaluated by proteomic and computational tools. Herein, we present a toolbox of fragments containing a 3,5-bis(trifluoromethyl)phenyl core that was equipped with chemically diverse electrophilic warheads showing a range of reactivities. We characterized the library members for their reactivity, aqueous stability and specificity for nucleophilic amino acids. By screening this library against a set of enzymes amenable for covalent inhibition, we showed that this approach experimentally characterized the accessibility and reactivity of targeted cysteines. Interesting covalent fragment hits were obtained for all investigated cysteine-containing enzymes.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery.

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Off-pathway 3D-structure provides protection against spontaneous Asn/Asp isomerization: shielding proteins Achilles heel.

Spontaneous deamidation prompted backbone isomerization of Asn/Asp residues resulting in - most cases - the insertion of an extra methylene group into the backbone poses a threat to the structural integrity of proteins. Here we present a systematical analysis of how temperature, pH, presence of charged residues, but most importantly backbone conformation and dynamics affect isomerization rates as determined by nuclear magnetic resonance in the case of designed peptide-models. We demonstrate that restricted mobility (such as being part of a secondary structural element) may safeguard against isomerization, but this protective factor is most effective in the case of off-pathway folds which can slow the reaction by several magnitudes compared to their on-pathway counterparts. We show that the geometric descriptors of the initial nucleophilic attack of the isomerization can be used to classify local conformation and contribute to the design of stable protein drugs, antibodies or the assessment of the severity of mutations.At any ­Asn/AspGly­ sites in proteins a spontaneous backbone isomerization occurs within days under physiological conditions leading to various forms of proteopathy. This unwanted transformation especially harmful to long-lived proteins (e.g. hemoglobin and crystallins), can be slowed down, though never stopped, by a rigid three-dimensional protein fold, if it can delay in the conformational maze, on-pathway intermediates from occurring.

The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein.

Invited for the cover of this issue is the group of András Perczel at Eötvös Loránd University, Budapest, Hungary and colleagues from Osaka University, Japan. The image depicts the amyloid buildup of an Exenatide derivate miniprotein (E5) monitored on a simplified hyperspace. Read the full text of the article at 10.1002/chem.201903826.

Structural impact of GTP binding on downstream KRAS signaling.

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states - an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg2+ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

1H, 15N backbone assignment and comparative analysis of the wild type and G12C, G12D, G12V mutants of K-Ras bound to GDP at physiological pH.

K-Ras protein is a membrane-bound small GTPase acting as a molecular switch. It plays a key role in many signal transduction pathways regulating cell proliferation, differentiation, survival, etc. It alternates between its GTP-bound active and the GDP-bound inactive conformers regulated by guanine nucleotide exchange factors and GTPase activating proteins. Its most frequent oncogenic mutants are G12C, G12D, and G12V that have impaired GTPase activity, thus induce malignant tumors. Here we report the resonance assignment of the backbone 1H and 15N nuclei of K-Ras wildtype, G12C, G12D and G12V proteins' catalytic G domain (1-169 residues) in GDP-bound state, and 13C of backbone and side chains of G12C mutant at physiological pH 7.4. Triple resonance data were used to get secondary structure information and backbone dynamics of G12C, the best-known drug target among K-Ras mutants. Simultaneous investigation of G12C, G12D and G12V mutants, along with the wild type form at the very same conditions allowed us to perform a comprehensive analysis based on the combined chemical shifts to reveal the effect of mutation at G12 position on structure. Intriguingly, the G12C and G12V mutants found to be structurally very similar at the three most important regions of K-Ras (P-loop, Switch-I, Switch-II), while the G12D mutant significantly differs at P-loop and Switch-II from the wildtype as well as G12C and G12V mutants. However, in Switch-I it hardly deviates from the wildtype protein.